Robust dynamical decoupling with concatenated continuous driving 
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The loss of coherence is one of the main obstacles for the implementation of quantum information 
processing. The efficiency of dynamical decoupling schemes, which have been introduced to address 
this problem, is limited itself by the fluctuations in the driving fields which will themselves introduce 
noise. We address this challenge by introducing the concept of concatenated continuous dynamical 
decoupling, which can overcome not only external noise but also fluctuations in driving fields that 
implement the decoupling sequences and thus holds the potential for achieving relaxation limited 
coherence times. We implement the scheme experimentally with nitrogen-vacancy (NV) centers in 
diamond, and demonstrate an improvement of the coherence time by an order of magnitude. The 
proposed scheme can be applied to a wide variety of other physical systems including, trapped atoms 
and ions, quantum dots, and may be combined with other quantum technologies challenges such as 
quantum sensing and quantum information processing. 



Introduction. — Coherent control of quantum systems 
has opened a promising route towards novel quantum 
devices for quantum technologies, such as quantum in- 
formation processing, quantum metrology and quantum 
sensing [IHl]. The performance of such quantum devices 
critically depends on coherence times of their constituent 
quantum systems which, in turn, is limited by uncon- 
trolled interactions with their surrounding environment. 
This results in a challenging but fundamentally impor- 
tant task in current quantum experiments, namely how 
to protect individual quantum states from decoherence 
by their environment while retaining the ability to con- 
trol the quantum dynamics of the system, in particular in 
solid state systems with characteristic complex environ- 
ments. Dynamical decoupling originated in liquid state 
NMR and allows for the resolution of complex spectra 
0- Significant progress has been made with the theo- 
retical proposals of various dynamical decoupling [6UT3] , 
and their experimental demonstrations p!4H23] . 

The recently developed dynamical decoupling schemes 
that require only continuous oscillatory driving fields 
[71 [5 , inherit the advantages of standard dynamical de- 
coupling, namely requiring no encoding overhead, no 
quantum measurements, and no feedback controls. More- 
over, they are easier to realize experimentally and are 
more naturally combined with other quantum informa- 
tion tasks, such as the implementation of high fidelity 
quantum gates pTHl4]. In principle, one can apply con- 
tinuous driving to reduce the noise suffered by a qubit 
considerably simply by increasing their intensity. Ran- 
dom and systematic fluctuations which are inevitably 
present in the driving field itself will ultimately limit the 
efficiency of dynamical decoupling. The deleterious effect 
of driving field fluctuations in particular will become sig- 
nificant when employing strong driving fields to achieve 
ultralong coherence times. Overcoming the limitations 
imposed by driving field fluctuations and thus extending 



coherence times further represents a key step towards 
the construction of quantum memory [24], highly sensi- 
tive nano-scale magnetometers [J-^ and error-resistant 
quantum operation [13 . It will also be of particular inter- 
est for Ti limited Rabi-type magnetic resonance imaging, 
the resolution of which highly depends on the stability of 
microwave driving fields [25 , 26 . 

In this Letter, we address this challenge by introduc- 
ing the scheme of concatenated continuous decoupling 
(CCD), which can significantly extend coherence times 
by protecting against driving field fluctuations. As a 
proof-of-principle, we implement a second-order CCD 
with a single NV center in diamond where a second, 
weaker driving field, reduces the impact of the ampli- 
tude fluctuations of the first order field. We demonstrate 
experimentally that CCD schemes with only a weak sec- 
ond driving field can already increase coherence times by 
an order of magnitude as compared to standard schemes 
based on a single drive. 

Concatenated continuous dynamical decoupling — We 
start by considering a two-level quantum system with the 
eigenstates |t) , ||), see Fig. 1(a). Its system Hamiltonian 
is Ho = ^ (It) (tl - \l) (II). Environmental noise causes 
fluctuations to the energies cj^, and thus the loss of 
coherence. To counter such effects, we can apply a driv- 
ing field on resonance with the energy gap between 
It) and It) as 

Hd^ = hQi cos {ujt)ax (1) 

where cr^ = It) (i| + |i) (t|- the interaction picture 
with respect to Hq and with rotating wave approxi- 
mation, we find H^^^ = ^ (1+) (+1 - I-) (-1), where 
|±) = ^ (It) ± It)) are the dressed states. In this ba- 
sis, the effect of dephasing noise now induces transitions 
among these dressed states, which are suppressed by an 
energy penalty as long as the noise power spectrum at the 
resonance frequency is negligible [14 . The decoupling ef- 
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FIG. 1. (Color online) a: Second-order concatenated con- 
tinuous dynamical decoupling: The first-order driving field 
with the frequency uj and the amplitude Qi creates the first- 
order dressed states |+), |— ), which suffer less from the de- 
phasing effect of environment noise, however are subject to 
the fluctuation in the amplitude Qi of the driving field. The 
subsequent second-order driving field with the amplitude Q.2 
(which is generally smaller than Q.i) is off-resonant with the 
detuning Qi and a relative phase cj) with respect to the first- 
order driving field. With such a second-order driving field, it 
is possible to overcome the intensity fluctuation of the first- 
order driving, and leads to the second-order dressed state W^) 
and \]^^) that are subject to the much reduced energy fluctu- 
ation, b: The diagram for the energy levels of the NV center 
electron spin. The NV spin triplet electronic ground state is 
splitted by an applied magnetic field. The effective two-level 
system used in our experiment is formed by the spin sublevels 
rris — (labeled as |^)) and rris — —1 (labeled as c: The 
procedure for the Ramsey experiment with the second-order 
dressed qubit: the NV electron spin is first polarized into the 
rris — sublevel, which is a coherent superposition of |tc/p) and 
l^c/p); the first- and second-order driving fields are simultane- 
ously switched on for time r followed by the optically readout 
of the population of = via spin-dependent fluorescence. 



ficiency will be limited if the noise has a wide range of 
frequencies, while it can be very efficient for slow baths 
e.g. in diamond [27 . The above analysis is based on 
the assumption that the amplitude of the driving fields 
is stable. In realistic experiments however, the intensity 
of the driving fields will fluctuate owing to limited sta- 
bility of microwave sources and amplifiers (the frequency 
instead can be relatively much more stable), and thus 
cause fluctuations of the energies of the dressed states. 
Achievable coherence times of the dressed qubit using a 
single drive are thus ultimately limited by the stability 
of the driving fields |25j , which appears as the fast decay 
of Rabi oscillation. 

The principal idea of CCD is to provide a concatenated 
set of continuous driving fields with decreasing intensities 
(and thus smaller absolute value of fluctuation) such that 
each new driving field protects against the fluctuations 
of the driving field at the preceding level. For example, 
fluctuations in the amplitude of the first-order driving 



can be suppressed by applying a second-order driving 
field as follows 

Hd^ = hrt2 cos [(cj + l^i) t + Lp\(Jx. (2) 

This second-order driving field is on resonance with the 
the energy gap of the first-order dressed state, see Fig.l 
(a), which actually describes rotation about the axis 
n((/p,— 7r/2) in the interaction picture (see SI for de- 
tails) and thus plays the role of decoupling the first-order 
dressed states from the fluctuation of Vti. With these 
two driving fields as Hd^ and Hd^^ we find the effective 
Hamiltonian in the second-order interaction picture as 

Hf^ = ^{\U{U-\U){U), (3) 

where the second-order dressed states are 

IU = /f (sin||t)+icos||;)) (4) 

IU = /|(cos||t>-isin||;)). (5) 

Therefore, one can encode quantum information in the 
above second-order dressed states as in Eq.(4-5). In gen- 
eral, we can apply n-order continuous driving fields on 
the condition that two subsequent drivings describe ro- 
tations about orthogonal axes in the corresponding inter- 
action picture. We remark that the relative phase does 
not need to take a specific value, once it is fixed the corre- 
sponding second-order dressed state can serve as a robust 
encoded qubit see SI for details. 

The decoherence of such n-th order dressed qubit stems 
dominantly from the fluctuation of the nth-order driving 
field and the residual effect of the fluctuations of the pre- 
ceding driving fields. As long as the orthogonality of two 
consecutive driving field is satisfied, the next consecu- 
tive field reduces the noise of the previous. Hence coher- 
ence times of the nth-order dressed qubit will mainly be 
limited by the fluctuations of the nth-order driving field. 
With a concatenated scheme in which subsequent driving 
fields have decreasing intensities, the effective dephasing 
will be sequentially suppressed and coherence times of the 
nth-order dressed qubit can be extended significantly. 

Experimental demonstration of CCD. — We have used 
NV centers in diamond to demonstrate the working prin- 
ciple and efficiency of our concatenated continuous dy- 
namical decoupHng scheme. As a promising candidate 
physical system for modern quantum technologies, NV 
centers have been used to demonstrate basic quantum 
information processing protocols [28^, as well as ultra- 
sensitive magnetometry and nano-scale imaging at room 
temperature j29H3T] . The electronic ground state of NV 
center is a spin triplet with three sublevels with mag- 
netic quantum numbers = and = ±1, and 
the zero field splitting is - 2.87GHz |321, see Fig.l (b). 
We apply an additional magnetic field along the axis of 
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(a) Single driving field 



(b) Second-order driving scheme 
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FIG. 2. (Color online) a: Coherent driven oscillation of NV 
center with a single microwave field. The Rabi frequency is 
0.1 — 40MHz. The blue curve is the experimental data, and 
the purple one is the fitting with a Gaussian decay envelope. 
The decay time is estimated to be about 2.5/is. b: Persistent 
Rabi oscillation by adding a second-order driving field, the 
intensity of which is 10 times weaker than the first driving 
(Qi = 20MHz). The four panels shows the very slow decay of 
Rabi oscillation beyond 1000/is. 



the NV center to split the energy levels of = ±1. 
The two electronic transition frequencies corresponding 
to rris = ^ rris = ±1 are determined to be 2042MHz 
and 3696MHz via an ODMR measurement. The effec- 
tive two-level system we use is formed by the sublevel 
rris = —1 = It) and rris — = Interaction with 
-^^C nuclear spin bath is the dominating source of deco- 
herence for ultrapure Ha type crystals which were used 
in our experiments. This relatively slow bath is ideal 
for decoupling experiments [21] . Significant progress was 
achieved in material engineering [33] but not all sources 
of magnetic noise can be eliminated in this way, in par- 
ticular for implanted NV defects [3T] . 

In our experiment, we first polarize the NV center elec- 
tron spin into the sublevel = with a green laser (532 
nm). We note that this is equivalent to the preparation 
of an initial coherent superposition for the relative phase 

7^ nTT, namely \ms = 0) = ||) ex cos ^ \\^) — sin ^ ||(^). 
For comparison, we start by applying a single driving field 
on resonance with the electronic transition rris = ^ 
rris = —1- We measure the oscillation of the population 
of state rris = 0- It can been seen from Figj2|a) that the 
Rabi oscillation decays in a few microseconds. By fitting 
the experimental data with a Gaussian decay envelope 
Si{r) = exp (— 6^r^/2) (see SI for details), we estimate 
that the decay rate of Rabi oscillation is bi ~ 88kHz. We 
stress that the fast decay is mainly due to intensity fluc- 
tuations of the microwave field, as the Rabi frequency is 
40MHz, which is strong enough to suppress the effect of 



magnetic noise almost completely (the nuclear spin bath 
formed by ^^C nuclei is slow compared to driving fre- 
quency). This statement will be supported and verified 
by the experimental data that we present here. 

To demonstrate the working principle of our concate- 
nated continuous dynamical decoupling, we add a second- 
order driving field with a weaker amplitude than the first 
driving. The measured population of the sublevel rris = 
can be written as (see SI for a derivation) 



^o(^) = ^ + ^ cos cos (1^1 r + (/:?) 

+ - sin (p sin {Qir + (p) cos (~^^ 



(6) 



We remark that it is not necessary to fix the relative 
phase (f between two driving fields for the estimation 
of coherence times. The measured signal averaging over 
random phase (p leads to Po{t) = § + | cos ((^ir) + 
I cos ((^it) cos (^r) . The oscillatory components and 
their decay features remain the same as a fixed phase. 

The second term in Eq.(6) represents Rabi oscillations 
from the first-order driving. Its decay is mainly due to 
the fluctuation of the first-order driving field, the effect 
of which now in the interaction picture is causing the 
relaxation of the dressed qubit and will be suppressed 
by the energy gap induced by the second-order driving. 
To demonstrate this, we first show that CCD scheme can 
sustain Rabi oscillations by applying a second-order driv- 
ing field, the intensity of which is ten times weaker than 
the first driving field. In Fig. 2 (b), we observe coherent 
Rabi oscillation after 300, 500, 800, lOOO/is, and find that 
it decays significantly slower than the one using only a 
single driving field (see Fig. 2a for comparison). Our ex- 
perimental data thus demonstrates that the effect of the 
fluctuation of the first driving field can be significantly 
suppressed by a second-order driving field. 
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FIG. 3. (Color online) a: The measured state population of 
the sublevel ms = as a function of the evolution time r. The 
blue curve is the experimental data, and the purple curve is 
numerical fitting. The Rabi frequency of the first driving is 
Qi — 40MHz, and the intensity of the second driving is three 
orders of magnitude weaker, b: The estimated coherence 
time of the second-order dressed qubit as a function of the 
amplitude of the second-order driving. The Rabi frequency 
of the first driving is = 40MHz. We fit the experimental 
data with a Gaussian decay envelop for the Rabi oscillation 
and the beatings from the second-order driving respectively. 
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We further estimate the dephasing time of the second- 
order dressed qubit by numericahy fitting the decay of 
beatings (i.e. the third term in Eq.(6)), which reflects 
the coherence of the second-order dressed qubit. The de- 
cay stems from the fluctuations of both the flrst and sec- 
ond driving flelds. Our experimental data, in combina- 
tion with numerical fltting, as plotted in Fig. 3 (b), shows 
that there is an optimal amplitude for the second-order 
driving fleld in order to achieve the best coherence time. 
By increasing the amplitude of the second-order driving 
one can suppress the effect of the fluctuation in the flrst 
driving fleld, however this will result in larger fluctua- 
tions of the second-order driving itself. The compromise 
between these two effects leads to the optimal choice of a 
second-order driving. In our experiment, with a second- 
order driving fleld about 1000 times less intense than 
the flrst-order driving, one can prolong the coherence by 
an order of magnitude from T2 = 2.5/i5 to T2 = 25/is, 
see Fig. 3 (a-b). Such a weak additional driving would 
not have apparent effect in suppressing environmental 
magnetic noise, nevertheless it signiflcantly improves the 
dynamical decoupling efficiency. This fact supports our 
claim that the fast decay of Rabi oscillation in Fig. 2 (a) is 
mainly due to the fluctuation in the driving fleld, which is 
shown to be a severe problem in conventional continuous 
dynamical decoupling with a single drive. 

We also demonstrate another important ingredient for 
concatenated continuous dynamical decoupling, namely 
flxing the relative phase between consecutive driving 
flelds, in order to generate a long living dressed qubit. We 
achieved this goal by using arbitrary waveform generator 
(AWG) to generate the flrst- and second-order driving 
flelds. In Fig. 4, we show an example measurement with 
an estimated flxed phase ip = 27r/5. 

Potential applications of CCD. — The CCD scheme can 
be combined with various quantum information process- 
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FIG. 4. (Color online) The measured state population of the 
sub level wis = as a function of the evolution time r. The 
blue curve is the experimental data, and the red curve is from 
numerical simulation with the estimated relative phase Lp — 
27r/5. The Rabi frequency of the first driving is Qi = 30MHz, 
and the second driving is Q2 = Oi/lO. 



ing tasks and make them robust against not only noise 
but also the fluctuations of decoupling flelds. The ini- 
tialization, coherent manipulation and readout for the 
second-order dressed qubit can be easily achieved with a 
continuous driving fleld at the frequency uOc = + (^2/2 
and with the same phase as the flrst-order driving. Such 
a qubit encoded with the sequential dressed states has an 
ultralong coherence time, and thereby can be exploited 
to construct a single-spin magnetometer to probe a weak 
oscillating magnetic fleld h(t) = Bag ^^^[{^i + ^/2)t]. If 
Bag is much larger than the fluctuation of the second 
driving fleld, the sensitivity will be ~ 1/a/7Y with Ti 
the coherence time demonstrated in Fig. 2(a); otherwise 
it scales as I/VT2 with T2 the prolonged dephasing time. 
The scheme can also be beneflcial for the construction 
of a precise noise spectrometer [26l |35]. We note that 
the inhomogeneity over an ensemble of quantum systems 
(namely spatial fluctuation) also leads to dephasing of 
the ensemble collective state. Our schemes can suppress 
this kind of spatial fluctuation and in the mean time can 
be robust against the inhomogeneity of the driving fleld 
itself, which may be useful in the ensemble-based mag- 
netometer and quantum memory [36 . 

Summary. — We have introduced the concept of CCD 
and implemented it experimentally. Using an NV center 
in diamond we demonstrated the superior performance 
of concatenated continuous dynamical decoupling com- 
pared to single driving flelds in extending coherence times 
of a dressed qubit. A qubit encoded in the concatenated 
dressed states is robust against both environmental de- 
phasing noise and intensity fluctuation of driving flelds. 
Our schemes can be applied to a wide variety of quantum 
systems where they may flnd applications in the construc- 
tion of nano-scale magnetometry and imaging e.g. with 
the NV center in diamond, and in the construction of 
fault-tolerant quantum gates that are protected against 
noise and control errors. 
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SUPPLEMENTARY MATERIAL 



Derivation of a second-order dressed qubit. — In our concatenated scheme, in general, we can apply n-order con- 
tinuous driving fields on the condition that two subsequent drivings describe rotations about orthogonal axes in 
the interaction picture. To illustrate the basic principle, we present the calculation details for the derivation of a 
second-order qubit. The system Hamiltonian is 

^^o = ^(it>(ti-i;)ai) (7) 

We apply two driving fields as 

= hVti cos {ujt)(Jx (8) 

Hd^ = MI2 cos [{uj + Vti) t + (^](Jx- (9) 

where = |t) (i| + |i) (t|- Iii the first-order interaction picture with respect to Hq and taking account of a rotating 
wave approximation, the effective Hamiltonian is 



(1) nrti 



It)ai+e^("^*+^)|;)(t|]. (10) 



The first driving field B.^^ describes rotation about the x axis. In the second-order interaction picture with respective 
to = ^^(Jx, the effective Hamiltonian becomes 

(2) ( . X .^^x 

Hj= (sm^^cTy — co^Lpdz) (11) 

where Gy and are Pauli operators in the spin up and down basis. Thus, one can see that the second driving field, 
although it is off-resonant with the original transition frequency cj, describes rotation along the axis n (9?, — 7r/2) that 
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is orthogonal to x. For each specific value of the relative phase (/?, the eigenstates of H\ ^ as in Eq.(4-5) in the main 
text can be used a dressed qubit to encode quantum information. In the above derivations, we have adopted the 
rotating wave approximations which hold when 1^2/2 <C l^i <C 2uj. These conditions are fulfilled in the present context. 

Rabi oscillation under a second- order continuous driving. — In our experiment, we initially prepare the NV spin in 
the state \ms = 0) = |^) oc cos ^ — sin ^ |^(^).The state evolution in the second-order interaction picture is 



cos2f/o(r) + i[l-/o(T)] 

eXp(ii^T)/i(T)/2(T) 



• 2 2 ^ 

- sm ^ COS ^ 



- sm ^ COS 2 



sm 



2 (p . 



2£exp(-i^T)/i(T)/2(T) 



/o(r) + i[l-/o(r)] 



(12) 



where the density matrix is written in the basis of {|t(^) ? lic^)}? foi'^) describes the relaxation from both the limited 
Ti of NV center itself and the fiuctuation of the first driving field, /i(r) and /2(t) represent the dephasing from the 
relaxation and the fiuctuation of the second driving field respectively. We remark that the first driving field is strong 
enough to eliminate the effect of magnetic noise in our experiment. We measure the population of the state \ms = 0) 
at time r, which can be written as 



Po(r) = Tr 

_ 1 

~ 2 " 



exp 



.^IT \ (2). . f.^lT 

-^^^^ jPi iV exp h^^a 



\i) ai 



^cos(pcos{niT ^ ip)fo{r) + ^sm(psm{niT ^ (f) cos (^r )/i(r)/2(r). 



In the above derivation, we have used 



exp ( ^^(^x I II) = cos + -) |t^) - sm 



(13) 
(14) 

(15) 



One can estimate the decay functions /^(t) without fixing the relative phase (p. The observable that we measure is 
thus the average of Po{t) over (p, which is as follows 



Poir) = \ + ^cos(i^iT)/o(T) + ^sin(i^iT)cos (^r \h{r)h{r). 



(16) 



The oscillation contrast will be half of the maximum contrast obtained by fixing a certain relative phase. If the second 
driving field is sufficiently strong, the effect of the first driving field fluctuation is suppressed and the relaxation 
decay /o(t) and /i(r) will mainly originate from the relaxation of the NV center itself. In such a case, even the oscil- 
lation of the third term of Eq.(8) may decay fast, we can still observe that the Rabi oscillations sustain for a long time. 



Numerical simulations of persistent Rabi oscillation. — We numerically simulate the exact dynamics with two driving 
fields, namely under the following Hamiltonian 



H 



-\- h [Qi + Si{t)] COS {ujt)ax + h [Q2 + ^2(^)1 cos [{uj + Qi) t + (p]o-x. 



(17) 



where = (|t) (t| ~ \i) ^i(^) ^^{t) represent the fiuct nations of the first and second driving fields respectively. 

We model the fiuctuations of driving fields 5i{t) i = 1,2 as an Ornstein-Uhlenbeck process with the correlation 
function Ci{t) =< Si{0)Si{t)) >= 6f exp( — |t|/rc.), where hi is the fiuctuation amplitude and Tc- is the correlation 
time of fiuctuation. We generate 2000 realizations of an Ornstein-Uhlenbeck process with the exact simulation 
algorithm [1 . We first numerically simulated the experimental data in Fig. 1(a) for the case of Rabi oscillation 
with one single driving field. The numerical result agrees well with the experimental data and the Gaussian fit 
of the decay envelope Si{r) = exp(— 6^r^/2) (see Fig. 2a of the main text). This supports that an Ornstein- 
Uhlenbeck process serves as a good model for the driving field fiuctuation. We estimate that the relative amplitude 
fluctuation of the driving fleld is about 2 x 10~^. We then use these estimated parameters for the numerical 
simulation of a second-order driving scheme. Fig. 2(b) shows Rabi oscillation with a second-order driving fleld 
using the same parameters as in Fig. 2(b) of the main text. We note that in the numerical simulation we flxed the 
relative phase, thus the contrast of Rabi oscillation is twice of Fig. 2(b). Nevertheless, the essential feature of the re- 
sults from our numerical simulation, namely the slow decay of Rabi oscillation, agrees well with the experimental data. 
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(a) Single driving field 
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(b) Second-order driving scheme 
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FIG. 5. (Color online) a: Numerical simulation of coherent driven oscillation of NV center with a single microwave field (red 
curve). The blue curve is the experimental data. The Rabi frequency is Qi = 40MHz, and its relative amplitude fluctuation 
is estimated to be (^i/Qi ^ 2 x 10~^. b: Numerical simulation of persistent Rabi oscillation by adding a second-order driving 
field, the intensity of which is 10 times weaker than the first driving (Qi = 20MHz). The plot shows that Rabi oscillation can 
sustain for much longer time than one single driving field. The four panels show the very slow decay of the Rabi oscillation 
beyond 300/is. The decay is mainly due to the relaxation of NV center itself, we use the value of the relaxation time Ti = 1.5ms 
in the numerical simulations. The relative phase is fixed as = 0. 

Coherent manipulation of a dressed qubit in sequential CCD scheme. — Dressed qubits can be coherently manip- 
ulated and coupled with each other via electron spin dipole interaction. To demonstrate this, we take = as an 
example. On can apply a continuous field with the frequency Uc = a; + 1^2/2 as 



He = ^cos {oJct)crx 



(18) 



In the first-order interaction picture with respect to the original Hamiltonian Hq = (|t) (t| ~ li) (ID? coupling 
can be written as 



hW = I 

2 



expH^t) It) ai+expH^t) i;) (tl 



In the second-order interaction picture with respect to i^Q^^ = ^^cFx^ it becomes 

rr(2) 9 f^'2.\ 
= :^ cos -—t ]cFx 
2 \2 J 

which induces Rabi oscillation between the second-order dressed states. 



(19) 



(20) 



The dressed qubits can be coupled via dipole-dipole interaction, e.g. between NV electron spins, as follows 

Hd-d = J{^SlSl — Si ' S2). 
In the first-order interaction picture, the electron-electron interaction in the driven Hamiltonian becomes [2] 



H 



(1) 



2 J S^S^ , 



(21) 



(22) 



The effective interaction between two second-order dressed qubits from the dipole-dipole interaction can written as 

J 



(2) 
d-d 



y^y 



1 2 



(23) 



where (Ty and are Pauli operators in the spin up and down basis. The effective entangling coupling between dressed 
qubits is thus feasible, and can be exploited to construct two-qubit gates. 
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